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TECHNICAL OBJECTIVES 

I. Quantum Calculations of Fluoroether Polymers 

The overall purpose of this research was to help understand polyfluoro- 

ether chemistry with a combination of experimental results and quantum 

mechanical calculations. Thus, we have worked in close colabor ation with 

Dr. Rosser at NASA-Ames Research Center who " contributed to the experimental 

counterpart of this research. The interpretation of the experimental results 

obtained by Rosser and coworkers about the relative thermal and hydrolytic 

stability of perfluoroalkylether substituted ljg^-oxadiazole 1 * 2 and 1,3,5- 

triazine 2 model compounds was the initial motivation of this investigation. 

The two general systems that we have considered are shown in figure 1. 

The perfluoro and perfluoroether substituents used in the experimental model 

compounds contain as many as seven carbon atoms. For the case of 1,2,4- 

oxadiazole, a model compound where R^ was a phenyl group bridging two 1,2,4- 

oxadiazole moieties has been synthesized and it has been found that its 

thermal stability decreases with respect to completed perfluorinated derivatives. 

it 

For the case of 1,3,5-triazine/was found that a branch (namely, a CHg group) 
at the carbon attached to the ring produced an extremely stable compound at 
high temperatures in different environments. Moreover, an unbranched 
perfluoro substituent was complete}^ hydrolized by water at 325C. 


Figure 1 
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(a) 1,2,4-oxadiazole derivatives 


(b) 1,3,5-triazine derivatives 



The theoretical aspect of this Investigation was mainly a study of the 
effect that different substitution produces on the electronic structure of 
the heterocycle. Changes in geometry, atomic charge populations, and molecular 
orbital energies were extensively examined. Substituents of smaller size 
than those present in the experimental model compounds have been used in the 
calculations. These substituents simulated the general, characteristics of 
the real ones, and have the advantage of reduced computational effort. 

XI, Quantum Chemical Calculations for a One-dimensional Conductor 

Another application of state-of-the-art quantum chemical methods under 
this contract was to a class of one-dimensional conductors. This class of 
conductors, of which. K^P^OT^BrQ^.xHgO (KCP) is one, has a chain structure 
consisting of staggered square planar tetracyanoplatinate complexes^. 

Different complexes can be formed with varying concentrations of 
halogens which modify the oxidation states of Pt from +2.0 to as high as +2.4. 
Examples of these complexes are given in Table X, A schematic diagram of 
the tetracyanoplatinate chain is shown in Pig. 2. For metallic conduction, 
the top valence band of these complexes must be partially occupied. Since 
the bonding along the chain is through d z * electrons, this implies that 
partial oxidation depletes electrons from the d z i band. Prediction of 

partial oxidation from specific orbitals requires an accurate knowledge of 

# 

the electronic structure. The objective of this research is to provide this 
knowledge of the electronic structure for the tetracyanoplatinate ion and to 
establish the orbitals from which partial oxidation occurs. The behavior 
of electrons in these orbitals are then correlated with the different Pt-Pt 
separations and electrical conductivities- 1 * 0 . A more fundamental understanding 
of the general trend of increase in conductivity resulting from increase 
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Table I 

Some Tetracyanoplatlnate Complexes * 


Complex 

Pt-Pfc (A) 

Conductivity -ft ^cm* 

Pt metal 

2.775 

105 

K 2 Pt(0M)4 

3.50 

10"7 

K2Pt(0H)4Bro.3 

2.89 

800 

0s 2 Pt(CH)4(IW) O-23 

2.872 

300 

C3 a Pfc(0M)4(FHP) 0i39 

2.833 

2000 

■» from Hefs (9) & (10) 




Figure 2 

[pt(CN)4 2 ] n chain 







in oxidation number for Pt coupled with decrease in Pt-Pt separation is sought • 


RESEARCH ACCOMPLISHED 

1, Quantum Calculations of Fluoroether Polymers 

As a computational tool we started using the semiempirical MINDO/3 
(Modified Intermediate Neglect of Differential Overlap) method of Dewar et alA 

(l) The s -triazine system 

MINDO/3 calculations on the unsubstituted 1,3,5-triazine yielded a planar 
molecule. This result is in agreement with experimental data**. Then we 
started studying the effect of monosubstitution of CF 3 CF 2 . The five fluorine 
atoms and the triazine ring were chosen to stay in a staggered conformation. 
Rotation of the triazine ring around the C-C bond was allowed, with the 
remaining geometrical parameters fixed. The goal of this calculation was to 

locate the neighborhood of the equilibrium conformation. The results of these 

* 

calculations indicated that the meet stable conformation was at a rotation 
angle of about 135 degrees. Consequently we started a calculation where we 
froze this angle and permit a partial relaxation of the geometry within the 
CF 3 CF 2 group. The results obtained look very unrealistic: a large deviation 
of the C-C-C and F-C-C angles from the tetrahedral value brings the fluorine 
atoms to an excessive proximity. 

Repetition of the calculation at 90° gave the same kind of anomalous 
result. We concluded that the cause of this problem resides in the under- 
estimation of lone pair - lone pair electron repulsion between adjacent 
fluorine atom, in MINDO/3. Subsequently we decided to use the different 
quantum mechanical method MNDO^ (Modified Neglect of Diatomic Overlap) which 
has proved to be more accurate than MINDO/3 for molecules containing hetero- 
atoms.^ >7 



I 


MNDO calculations were then performed for vmsubstltuted s-triazine, 
and' for mono substitution of CF3CF2 . For the latter case the optimized 
geometry was much more realistic, and the heat of formation obtained was 
much less exothermic than that obtained with MINDO/3. For unsubstituted 
s-triuzine, MINDO/3 and MNDO predicts practically the same geometry for the 
ring. This agreement would have been expected since the absence of fluorine 
atoms reduces the error coming from the underestimation of lone pair - lone 
pair electron repulsion. 

MNDO calculations with the substituents CF(CF3)CFg and CF(CF3)OCF3 were 
also performed. We successfully optimized all the geometrical parameters. 

No significant differences at the triazine ring are detected in geometries, 
charge distributions and molecular orbital diagrams for the three different 
substituents employed. 

A close scrutiny of the molecular orbitals indicated that the loss of 
basicity of the triazine nitrogen atoms upon substitution is due to a drop 
of 0.8eV in the energy of the highest occupied molecular orbital of 
symmetry®'. Nevertheless differences in basicity between the three substituted 
s-trlazine is not apparent. The differences in hydrolitic stability could 
be of thermodynamic and / or kinetic nature. Both possibilities are currently 
being examined. 



We started studying the 1,2,4-oxadiazole simultaneously with the s-triazine 
system. Consequently some MINDO/3 calculations were initially performed. 

MINDO/3 results for the unsubstituted heterocycle revealed that all the atoms 
of the ring lie in the same plane in agre»2ment with experimental data®. 

A OFg group was substituted at carbon 3 (see figure la), and rotation of it 
aroimd the C-C bond was allowed. The results indicated that: (a) a free 
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rotation of CFg is expected at roam temperature, and (b) the rotation is 
practically rigid. 

Due to the shortcomings of MINDO/3 as seen in the s-triazlne calculations 
we performed MNDO calculations on the unsubstltuted 1,2,4-oxadlazole in 
order to see how the results differed from those of MINDO/3. She MNDO heat 
of formation was about 17 kcal/mol above the corresponding MINDO/3 value, 
suggesting that some underestimation of electronic repulsion might be present 
in the MINDO/3 calculation. Also the geometry was somehow different. 
Nevertheless a planar molecule was predicted in agreement with MXNDO/3 result. 
With these results at hand, and recognizing the improvements of MNDO over 
MINDO/3 for system with heteroatome^*?, we decided to continue the study of 
this system with the MNDO method. 

We have performed MNDO calculations for a series of molecules with 
different substituents and R 2 (figure 1(a)). Table II shows all the 
different cases considered. No details of all the different calculations 
are presented here. For this purpose the draft of a manuscript that we are 
preparing for publication is attached as an appendix. Here, only a summary 
of the conclusions achieved is given. Two important factors determine the 
thermal stability of the 1,2,4-oxadlazole structure: (a) the total negative 
charge ipside the ring, and (b) the degree of homogeneity of such charge 
distribution. For phenylated derivatives, a non-planar conformation is 
the energetically most favored. Increasing the temperature would start 
populating the planar conformation, and as a result a greater number of 
electrons would be accumulated inside the oxadiazole ring, thus leading to 
reduced thermal stability. Stabilization of phenylated derivatives is 
suggested by avoiding the planar 1 conformation as well as by the introduction 
of electron-withdrawing groups at ortho and para benzene ring positions. 
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Table II 


Ri 


R 2 


H 

CHo 

CF 


c 2 F 5 

CS'3 

c 2 f 5 

c 6 h 5 * 
C6H^F 2 + 

cf 3 


H 

H 

H 

H 

C? 3 

CFo 

if 

CFo 

c 6% 



* C6H5 is 
+ c 6«3 f 2 is p u 

II. Quantum Chemical Calculations for a One-dimensional Conductor 

The proposed calculations are performed using the self-consistent- 
field-Xrt- Dirac-Scattered-Wave (SCF-Xo(-DSW) method"^* for the clusters 


[pt(0N)4]' 


r t 


-4 


and |pt(CN)4| 3 for several Pt-Pt separations. This method has 
proven its utility in treating systems containing transition metals, and it 


treats relativistic effects exactly, which is vital to an accurate description 
of the Pt-Pt bond. The density of states is also calculated for each case. 
This study is intended as the first of a series of electronic structure 
calculations for one-dimensional conductors. 

The monomer and trlmer clusters each has point group symmetry. The 
density of occupied states for the monomer jpt(CN)4j" 2 is shown in Fig. 3. 

The orbital with d z %character is in the middle of the Pt-ligand valence 
complex. This is to be expected since the monomer does not represent the 
extended chain system. 



* Density of occupied electronic states for PtCCN )^" 2 based on the 
• SCF-X^-DSH molecular orbital energies 
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Density of occupied electronic states for jptCCN)^ 6 calculated from 
the SCF-X0f-3>SW molecular orbital energies 
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Substantial code developments have been made for the SCF-X«T-DSW computer 
programs, so that It can handle the 27 -at cm trlmer cluster ^t(CN)lJ 3 . 

The calculation for a Pt-Pt separation of 3.1A has dust been completed and 
the density of states Is shown in Fig. 4. Preliminary analysis reveals that 
the antibonding d^ states are located at the top of the valence complex* 
However, the density of states near the valence band edge (Ey) is not 
exceedingly high. This is reasonable since the Pt-Pt separation used falls 
between those for the insulator and conductor compounds as seen in Table I. 

It is fair to speculate at this point that the density of states near Ev 
will be lower for large Pt-Pt separations (insulator) and much higher for 
smaller Pt-Pt separations (conductor). More detailed analysis of the trimer 
electronic structure as well as calculations for three other separations 
are ongoing. 
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MXNDO/3 and MNDO calculations on Polyperfluoroalkyl Compounds, J.P. Lopez, 
C*,Y. Yang and R, Rosser, The Third Annual West Coast Theoretical Chemistry 
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■ MnUrniir nrhifl i Ttininy A pplied to Structural Changes In 
1,2,4-Oxadiazole Substituents 


INTRODUCTION 

The literature describes the relative thermal stabilities of various 

1 2 

substituted 1,2,4-oxadiazole conpounds, ' • General organic chemistry theory 
leaves many questions unanswered concerning the rationale of existing experi- 
mental data. 

In this paper a theoretical study of the effect that the electron donor 
power of different substituents at positions 3 and 5 of the 1,2,4-oxadiazole 
ring lias on the electronic structure of the heterocycle is undertaken, The 
results reported are all from the semi-enpirical quantum mechanical Modified 

3 

Neglect of Diatomic Overlap (MNDO) method , This method has proven to be 

4 

successful in the description of organic molecules containing heteroatoms , We 
will discuss the results of Mullilcen population analysis and a correlation 
molecular orbital diagram for a variety of selected molecules , 

Figure I shows the general molecular system for which several calculations 
lrnve been performed and, for simplicity, it is designated R^-OXAH^ in this paper. 
Tne objective of the work was to understand how the different substituents and 
1*2 modify the thermal stability of the 1,2,4-oxadiazole ring. We seek an answer 
to the question of how a particular substituent affects the electronic structure 
of the heterocycle, that is, does it favor dissociation of the ring or stabilize 
it. 

Different combinations of the following substituents have been employed in 
the calculations: hydrogen, methyl jperfluoromethyl, perf luoroethyl , phenyl and 

. wbi -fluorophenyl . 

V ' cL.liuf}' 


/the 
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iThg.rhalra-ef the substituents h avo b een ■ m otei va feadJs^fcha ■ a wi« t e ne e - o f 
similar .experimental model compounds that are relevant to the synthesis of 
thermally stable polymeric materials, P ar t icul arly, Experiments have revealed 
that substitution of perfluoroalkyl groups on both carbon positions of the 
oxadiazole structure greatly stabilizes the ring system^ j' Jtitfhe introduction 
of a phenyl group on the oxadiazole ring yields c^npounc \a which are not as stable 

as the completely perfluorinated derivatives Vtfc r 1 f? vm ■pw' (74w» the'v^lporinated 

bM 


derivatives, it is known that the unsubstituted 1,2,4-oxadiazole is^veffy- unstable 

> . v’U v 

Ppiis stability is significantly enhanced by replacement of one of the hydrogens 

£ 

by a methyl group:. 

Changes in atomic charge populations, geometries, and molecular orbital 
energies were extensively examined and in particular the effect of rotational 
conformers on the thermal stability of the phenyl substituted oxadiazoles will be 
discussed. 

For ' calculations where the phenyl group was present, the geometry 

p 'vjas kept fixed at that of free benzene molecule. A complete 


optimization of the 




geometrical parameters washer formed. For the 


cases where the phenyl group ye absent , a complete optimization of the atomic 
arrangement was undertaken. 

RESULTS 

(a) Rotational Barriers 

The calculations on all th a diff a r a nte mo loculi R^-O^Ol-R^have shown 
that, the five atoms of the oxadiazole ring lie in practically the same plane. 

Rotations around single bonds present insignificant barriers, except for the 
case of rotation of the phenyl group. Thus, the rotational barrier for 
CF^-QXA-CF^ is only 0.1 keal/mol, whereas for C^K^-OXA-CF^, C^H^-OXA-H, and 
CF^-OXA-CgH^ , they are 1.34, 1.28, and 1.16 keal/mol respectively. 
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The potential energy curve for rotation of^phenyl group around the 


C-C single bond ip depicted in Figure XX for CJH,.-0XA-CF o . The two minima 
are at 64° and -64° away frcm the planar conformation. Since these t^LfiLd 
conformations are separated frcm the plana/ conformation by 1.34 kcal/mol, it 


1 ~ '■**•**»►: *’ \Wf\~ rp**"*w 

. is expected that at room temperature pjaafafertAy* nnly the tyiaWd conformations 

f'W «rU fojtfid? u 7)'rL 1 *"* . *»PUmr 

va*e populated*, A quasi-free pass betoken the two oJtofrd conformations must be 

A 

expected because of the small barrier of separation between them (only 0 . 2 kcal/mol) , 
(b) Mulliken Population Analysis 

Results within the 1,2,4-oxadiazole ring are presented in Table I. 

A* 

\ diffnr^nt way to visualize the effect of substitution on the rearrangement of 
the atomic charge , is by examining the electron donor power of the substituents 
at positions 3 and 5 of the 1,2,4-oxadiazole ring. These results are shown in 
• Table III. Inspection of the table reveals that the electron releasing power 
increases according to the following order: 

C 6 H 3 F 2 (NPL > < C 2 F 5 < GF 3 < C 6 H 3 F 2 < C fi H 5 (NFL ) 


vu* r ?Uvi<tv 


<C 6 H 5 (pfc) <H<CH 3 


VdWlT 


The general ordering is in agreement with could be expected from 
qualitative argunents. It is interesting to observe the effectiveness of 
fluorine atoms in reducing the electron releasing power of the benzene structure. 
Also, it is seen that position 3 is a better s.mlr nf plftntaens than position 5, 


l,*H\ jk 


probably due to the o ^ fc i a gt^pr'ia a uu^j> of the more electronegative oxygen atom 
(conpared with nitrogen) . 

(c) Correlation Molecular Orbital Diagram 

Diagram for seven of the molecules of Tables I and II is depicted in 
Figure III. In addition, and for reference purposes, the benzene molecule has 
also been included. The diagrams show the energies of only the highest double 
occupied molecular orbitals of each molecule. In principle, a molecular orbital 
is delocalized throughout the entire molecule; consequent^it should have contributions 
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from the oxadiazole ring as well as from the substituents. Nevertheless, an 

inspection of the highest occupied orbitals reveals that they are practically 

» 

restricted to the oxadiazole and/or the benzene ring for all but one of the cases 
examined (the exception is CH^-OXA-H) , Molecular orbitals with significant 
contributions of the perfluoroalkyl substituents are much lower in energy. 

Since the 1,2,4-oxadiazole ring lies practically in the same plane for all the 
calculated structures, it is possible to classify the molecular orbitals that 
are "localized inside the ring" as of a and n synmetry. This notation is used 
in Figure III for oxadiazole [IT (o) and cr (o) ] and for benzene [II (b) and a (b) ] . , 
This notation is an approximate representation of the ID because of the 1,2,4- 
oxadiazole ring yf no^pefeectly planar, and because some tails from the 
perfluoroalkyl substituents are present in these orbitals. Moreoever ^the notation 

is useful because it shows the gross general character of the molecular orbitals. 

* 

The MO diagrams of (NFL) -OXA-CF^ and (NPL) -QXA-H correspond 

to conformations where the oxadiazole is perpendicular to the benzene ring. The 


reasons why these conformations have been considered instead of the most stable 
taalk^^ifl^ue^due to: (i) an easier assignment of the character of the 

irolecuiar orbital is performed, and (ii) the energies of the highest occupied MO 
are practically identical in both kinds of conformations. 

(d) Geometry 

A complete optimization of the geometry within the 1,2,4-oxadiazole 
ring was performed in all cases . In general a very weak dependence of geometry 

upon substitution was found. As an illustration of die rearrangement of the 
atomic positions within tlie oxadiazole ring, the results of the eo/uilibrium 
distances bond angles for four selected molecules are shown in Tab.ie III . 
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DISCUSSION 


'-VV 


(A) Perf] varoalkyl Derivatives 

Degradation studies have shown that 3 , 5 -bis (perf luorog^'i&ptyl) -1 ,2,4- 

I 2 

oxadiazole is thermally stable ' . It remains undissociated for temperatures 

as high as 325°C in nitrogen and oxygen atmospheres ; however, replacement of the 
perfluoro-jj^heptyl substituent at position 3 by a phenyl group yields a structure 
which dissociates completely at 325°C, 

Comparisons of the calculations for CF^-OXA-CF^, ^^“OXA-CF^ and 
CgHjj-O^-CF^ help to rationalize the experimental results. For the phenylated 
derivatives there is an obvious relationship between the heating process and the 
distribution of conformational populations according to the results in Figure 
II. A rise in temperature increases the relative population of the p-OP ar 
conformation; whereas j ^ow~ tenperatur^ t most of the molecules will exist in the 
nonplanar conformation^" It is worth noting that at 325°C the thermal energy is 
of the same order of magnitude as the calculated rotational barriers between both 
extreme conformations; consequently^ a significant increase in the population 
of the planar conformation is expected at this tenperature. The immediate 
question that arises now is whether or not the existence of the planar conformation 
Soll^be responsible for the loss of stability of the 1 , 2 , 4-oxadiazole ring. Thus 
we need to see what properties of the heterocycle are modified in the transition 
from the non-planar to the planar structure, and then to see how these changes 
could eventually lead to a ring collapse. Changes in geometry, atomic charge 
population and molecular orbital diagram will be examined. 

(i) Geometry: Even when the presence of the benzene ring, particularly 

in the planar conformation, produces the largest relative is'in in the 

oxadiazole ring, the changes are so small that it is difficult to know if they are 
related to ring rupture, 
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(ii) Atomic Charge Population: As can be seen from Table II, the 

planar conformation yields a greater accumulation of electrons inside the 

* • 

oxadiazole ring. This excess of negative charge could produce strain inside the 
ring leading to an eventual dissociation. On the other hand, the stable conpounds 
CF^-OXA-CF^ and C 2 F^ -OXA-CF^ present the least accumulation of negative charge 

in the ring. In suimary, the relative thermal instability of the phenylated 

\ 

derivatives seems to be di£ to the greater accumulation of electrons inside the 
1,2,4-oxadiazole ring produced by population of the planar conformation at high 
temperatures . 


(iii) Molecular Orbital Diagram: The ionization potentials obtained , 
from Figure III, and assuming validity of Koopman'a theorem, are: (p^-OXA-CF^ 

(9.92 eV), C 6 H 5 (NPL)*OXA-CF 3 (9.99 eV), and CF 3 -OXA-CF 3 (12.9 eV) . The higher 
ionization potential of CF 3 -OXA-CF 3 could explain the relative high oxidative 
stability of the completely perfluorinated derivatives. Introduction o ^benzene 
ring decreases tlie ionization potential, so the highest occupied MO's have electrons 
that are more readily transferable to an oxidizing agent. Inspection of MO 


diagram (Figure III) shows that the two highest occupied MO for the 
conformation of C^-OXA-CJ^ are mainly localized in the benzene ring (at -9.99 
and -10.10 ev) . The next highest occupied MO belongs (pra ctica l!^ to the oxadiazole 
ring with an eigenvalue of -11.77 eV. For tlie planar conformation, as 
be expected, a delocalization of the electronic charges to both ring systems is 
found. Even when the IP does not change in a considerable amount, the character 
of tlie highest MO is altered, Now the KCM) (at -9,92 eV) has a significant contri- 
bution from the 2 p z orbitals of the atoms located at the oxadiazole ring. For 

the planar conformation tlie oxidation itJylome extent involves removal of an 

1 i A > 

electron from the oxadiazole. This would explain the decrease in oxidative 
thermal stabi3,ity of the phenylated derivatives: The probability of oxidation of 

the oxadiazole ring increases with an increment of the temperature because the 
planar conformation has a greater population. 
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(B) Establishment of a More General Relationship 

The number of electrons accumulated Inside the 1,2,4-oxadiazole ring 

t 

can be viewed as an "index" to predict differences in thermal stability. But, 


although this irjtex explains the differences in stabilities between the completely 

Vl ktV'l , . , , St- 

ranin the stabilities of ^ 


fluorinated and the phenylated derivatives , )\ 

^QyOXA-H and H-QXA-H aVe c onpated./ CH^-OXA-H is thermally more stable than 
H-OXA-H despite the fact that the oxadiazole ring has a greater negative charge 
for the former case (see Table II) . Then a better t ndewnoo d g t i n li e frwnd . ■fti nh 
}Mfi?^n^ex should include (in addition to the total number of electrons released 

to the oxadiazole system) a measure of the charge distribution around the ring. It 
could be intuitive that an homogeneous distribution of a definite number of electrons 
will produce less strain in the ring than the same number of electrons but heterogenously 
distributed. One way to have an indication of the homogeneity of the charge 
distribution is looking at the standard deviation of the charge, The results for 
the standard deviation and the new tentative index ft are presented in Table IV. 

The index is defined asX“ total number of electrons released to oxadiazole 

OL . 

ring X standard deviation. Thus irr smaller value of iL.is associated with a more 


stable compound. A perfect correlation between^ and those structures for which 


experimental information is available found. 

(C) Predictions of Thermal Stability 

The‘ index/l can be used to predict the relative thermal stability of 
compounds for which experimental information is not available. The different 
values of-Z^are shown in Table IV. 

The thermal stability of the phenylated derivatives can be improved 
by decreasing the electron donor power of the benzene structure. This goal can 
be accomplished by attaching electron withdrawing groups to the benzene ring, and 
by avoiding the planar conformation, Both effects are nicely seen in C^H^-OXA-CF^ 

.hi . 

which is predicted to be ^'rnore stable structure (the smallest./!) . In this case, 
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a veryJiig' rotational barrier of 5.8 kcal/mol separates the non-planar (now at 
78°'. from the planar conformation, A temperature as high as 2630°C is needed in 
order to provide a thermal energy of 5.8 kcal/mol, thus the non-planar conformation 
is expected to have greater population. 

Examination of M3 diagran for C^y^-OXA-C^ reveals a general lowering 
of the energy levels compared with C^-OXA-CF^ . The ionization potential 
increases from 9.9 to 10.3 eV, In addition, the character of some of the highest 
occupied MD’s is also modified. For the non-planar conformation the two highest 
occupied MD's are mainly of benzene character (at -10,27 and -10.80 eV); but the 
next highest MO (at -11,84 eV) has contributions from both rings in contrast to 
the equivalent MO of (NPIO-OXA-CF^ which is only localized in the oxadieole 

^ y J a V 

)f tUeSxaQiazyie iin^ rnSrief o; 
ivez Jit r^not e^ect^fc! oVat a siafificaijcpgpuJlatibn 

v / 

gre^oer ryertnes* v> regcva/ or ej^eotrons rrom zne 's^aiazoig' rju 

"* The value of the total charge released to the oxadiazole ring is 
about the same regardless if the- phenyl group is attached to carbon 3 or carbon 5 
as can be observed frcm Table TV for C^H^-OXA-CF^ anc * CF^-OXA-CgH^ respectively. 
ButJ*phenyl group at position 5 generates a less homogeneous distribution of this 
charge as can be deduced from the values of the standard deviation shown in Table 
IV. Accordingly, substitution of ^phenyl group at position 5 o /l,2 ,4-oxadiazole 



predicts conpounds thermally less stable, 


Consequently, a greater inertness to removal of electrons from the 
oxadiazole ring is expected for the non-planar conformation of C^H^-OXA-CFj. 

This relative inertness is partially lost for the planar conformation of 
C 6 H 3 F 2’ OXA " CF 3 k ecause a larger contribution of the oxadiazole ring in one of 
the highest occupied M)’s is observed. (This would be expected due to it delocali- 
zation to both rings - see Figure III) . Accordingly, more electrons of the 
oxadiazole ring are now available for removal. Nevertheless, due to the huge 
rotational barrier, population of the planar conformation will be important 
only for very . igh temperatures (over 2000°C) , In summary, the oxidative thermal 
stability of the oxadiazole ring is predicted to be improved by the presence of 
the 2 -fluorine substituents on the phenyl ring. 
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SLEMARY 

Hie loss of stability of the 1,2,4-oxadiazole ring with various substituents 
is due to population of the less stable planar conformation at high temperatures, 
which by accumulation of excessive nurtoer of electrons inside the heterocycle could 
lead to dissociation. The thermal oxidation seems to be favored by II delocaliza- 
tion since sane fraction of the electrons are now located in high energy levels 
of the 1,2,4-oxadiazole ring, 

Recognizing that two important factors determinate the stability of the 
1,2,4-oxadiazole structure* p* the total negative charge inside the ring, and its 
corresponding distribution, a tentative index to predict thermal stability has 
been proposed, 

Stabilization of the 1,2,4-oxadiazole ring for the cases of phenylated 
derivatives is suggested by avoiding the planar conformation as well as by the 
introduction of electron witdidrawing groups atT)ef?zene ring positions , 
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TABLE I. Mulliken Population Analyses Within 1,2,4-Qxadiazole Ring. The charge on each atom is given in 
muter of electrons. The mufcering of the atoms is that given in Figure I. 
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TAELE II. Nunber of electrons released to the 1,2,4-oxadiazole ring by several substituents . 

For the sake of clarity, die substituent group is show in the parenthesis next to 
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TABLE III. Equilibrium Geometry of 1,2. A-oxadiazole ring 
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TABLE IV. Standard deviation of the atomic charge distribution inside the 1 , 2 , 4 -oxadiazole ring 
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Bimonthly Technical Progress Report X 
for the period 

October 1, 1981 to November 30, 1981 


"Quantum Chemical Calculations for 
Polymers and Organic Compounds" 


I. Tetracyacoplotinate Compounds 

The SCF-X<rf-DSW relativistic calculations for Pt(CN)^^'’ that we completed 
previously were performed with the following maximum 2-value for partial-wave 
expansion at each center: Outer sphere, * 4; Pt, * 3; C and N, 

l max ■ 2, A non-relatlvistic calculation with an identical basis set was 
undertaken in order to establish a detailed comparison with the relativistic 
results. This calculation was performed also with the SCF-Xrf-DSW method, but 
setting the speed of light to infinity. The results should be identical to 
those obtained using the SCF-X<*-SW (Schroedlnger) method. However, the molecular 
orbital energies resulting from this non-relatlvistic calculation are signifi- 
cantly different from those obtained by Interrante and Messmer in a previous 
SCF-Xrf-SW calculation. 1 There are two fundamental differences between our 
calculation and that of Interrante and Messmer: (1) they used a smaller basis 
set (outer sphere, Jl max - 3; Pt, A jnax - 2; C and N, - 1); (2) in our 

calculation the orbitals 5s and 5p on Pt are considered as core states, whereas 
Interrante and Messmer Include these among the valence states. Consequently, 
we have a total of 48 valence electrons instead of 56. In order to have an 
indication of the magnitude of the ligand field shift and splitting of the 5s 
and 5p orbitals of Pt, we performed an energy search for them by using the 
last converged potential obtained for our non-relativistic calculation of 
PtCCN)^-. We include also in this search the 4f states since they are even 
higher in energy than the 5s orbital. In all cases the splittings and shifts 
were negligible. Thus, we feel that the consideration of 5s, 5p and 4f as core 
states would not produce any significant change in the SCF energy of the 
molecular orbitals. In order to understand the Importance of the size of the 
basis set we performed a non-relatlvistic (c-*cc) SCF-XjC-DSW calculation with 
identical Amax values than those used by Interrante and Messmer (however, we 
kept only 48 valence electrons). The molecular orbital energies are only 
slightly different from those obtained in the non-relativistic calculation with 
the extended basis. The difference with Messmer calculation is still significant. 
Thus, in order to ascertain the proper convergence of the DSW results to the 
Schroedinger SW ones, we performed a calculation with the latter. The results 
obtained are identical to those of DSW with c^&o. We think that the error in 
the Interrante and Messmer calculation could come from three different sources: 

(1) computer precision, (2) convergence criterion, (3) since overlapping 
interatomic spheres are considered, negative intersphere normalization might be 
present in their calculation. 
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Finally, a relativistic SCF-X«M)SW calculation with a basis set of 
intermediate size was also performed. In this calculation of outer sphere 

was 4, H max of Pt was 3 , and of C and N was 1. Molecular orbital energy 
diagrams for all the calculations already described are shown in figures 1 and 
2. The heading of each column has the following meaning: Measmer 56 corresponds 
to the Interrante and Mesamer calculation; Schro48 is a SCF-36^*SW (Schroedinger) 
calculation with the reduced basis set; NRDlrac48 is a non-relatlvistlc SCF-X«*-DSW 
calculation with the reduced basis set; NRDiPoFnAS is a non-relativistic 
SCF-Xtf-DSW calculation with the extended basis set; ReDiPoFn48 is a relativistic 
SCF-XuC~DSW calculation with the extended basis set; and 04Pt3CNl is a relativistic 
SCF-fcfc-DSW calculation with the basis set of intermediate size. 
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II, Fluoroether Compounds 

A MNDO calculation for the protonated specie shown in figure 3- has been 
completed, The convergence was very slow: complete optimization was obtained 
after 95 SCF calculations. The protonation appears to a highly endothermic 
process, with a protonation enthalpy of the order of 200 kcnl/mola, Addition 



CF 3 


of a proton produces significant changes in the geometry of the triazing ring; 
the bond distances as well as the bond angles are modified. The molecular 
orbital energies are all smaller compared with the unprotonated structure; a 
manifestation of this is the increase in the ionization potential from 12.3 to 
17.2 eV. 


The version of the MNDO program that we have been using until now can 
handle a maximum of 75 basis orbitals. We have been working in modifying the 
code with the aim to study systems of larger molecular size. Recently we have 
succeeded in this task: molecular systems with up to 160 basis orbitals (only 
valence electrons are considered) can now be considered. We have tested the 
program with the molecule shown in figure 4 which contains a total of 156 orbitals. 


Figure 4 


CFg — 0- 



We know positively that the iterative process is working correctly; we have been 
able to perform one calculation with 99 iterations; nevertheless self-consistency 
has not yet been achieved. The CPU time used in this calculation was about 
5 hours (in a VAX machine) . We have incremented the number of iterations but for 
some unknown reason the execution is sometimes interrupted or the output file is 
empty. Thus, clearly there is a problem with the computer system that we are 
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using. This machine is very slow, and when the elapsed time is too large (for 
a GPU time of 5 hours the elapsed time was 17 hours) there are inconsistencies 
in the response of the system. These facts are auggejting that MNDO calculations 
for large systems in a VAX computer are becoming to be very impractical. Thus, 
we are planning to transfer our modified MNDO program to the CRAY 1 computer 
which should be able to handle these systems more efficiently. 


References 
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Bimonthly Technical Progress Report IX 
for the period 

November 23, 1981 to January 22, 1982 

"Quantum Chemical Calculations for 
Polymers and Organic Compounds" 


I. Tetracyanoplatinate Compounds 


Both relatl v^tlc^ and^^ on the monomer have 


contributions front'd -like orbitals} 37% (d xy component) for non-relativistic 
and 37% for relativistic. In both cases there is an orbital that is mostly 
platinum d z 2-like. This orbital is 0.48 eV and 1,13 eV below the HOMO for the 
relativistic and non-relativistic calculations respectively. The d z 2 character 
also appears in some other orbitals for both kinds of calculations, 
la 




The oasis sec used m these two calculations corresponds to the one that we 
called of "intermediate size" in the previous report, that is JU 4 for the outer 
sphere, Jt$3 for Pt atoms, and JL$1 for C and N atoms. Energy level diagrams 
for these two distances are shown in figure 2. The relativistic energy levels 
of the monomer (obtained with an equivalent basis set) have also been added to 
figure 2 for comparison purposes. For the calculation at a Pt-Pt distance of 
3.5 A we found that the HOMO has 427, Pt character. Most of this charge is on 
the central Pt, having an Important contribution from d orbitals in the plane of 
the ligands. The orbital with the most d z 2 character (49% at the central atom 
and 57o at each peripheral Pt atoms) appears 1.28 eV below the HOMO. When the 
Pt-Pt distance is shortened to 2.89 A, the HOMO has now a greater amounc of d z 2 
charater: 57.5% (30.1% from central Pt and 13.7% from each peripheral atom). 

In addition, this HOMO is much more localized on Pt atoms than the HOMO at the 
larger distance: 857. (at 2.89 A) vs 427. (at 3.5 A). The d z 2 character is also 
present in some other orbitals for both distances. Its partial contribution 
to each MO can be easily seen in figure 3 where the d z 2 partial density of states 
of the central Pt atom is shown for each separation for a wide range of energies. 
These results so far are pointing towards the following preliminary scenario on 
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the relationship between partial oxidation and conduction mechanism, Tha 
oxidizing agant (tha halogan) la removing alatrona from planaa parallel to tha 
Pt-Pt chain; subsequent ahortening of tha Pt-Pt dlatanca producaa a rearrange- 
ment of the orbltala which reaulta In the partial occupation of a d x 2-lik« 
energy band. Thla partially filled band would be responsible for conduction 
along the Pt-Pt chain. We alao performed calculatlona where 0.9 electrona 
(In effect creating an oxidation atate of +2.3 for each Ft atom) were removed 
from the d z 2 like HOMO. In thla caaa the HOMO atlll conalata of mostly d s 2 
character, but the d 2 partial danalty of atate showa a algnlf leant decreaaa 
at the top of the valence band, Extenalva analyala of thaaa reaulta are ongoing 
and will be preeentad In the next progreae report. 

The reaulta obtained for the monomer and trlmer will be praaented at the 
March meeting of the American Phyalcal Society to be held In Dallaa, Texaa, 


II. Fluoroether Compounda 

A manuacrlpt entitled "MNDO Calculatlona Applied to Structural Changea In 
1,2,4-oxadlazole Subatituenta", la being Bent to Journal of Molecular Structure 
for publication. 


References 
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Bimonthly Technical Progress Report III 
for the period 

January 23, 1982 to March 22, 1982 

"Quantum Chemical Calculations for 
Polymers and Organic Compounds" 


The results of several calculations for the trimer (Pt(CN) . 2 ^ 3 have 
been extensively analyzed, and a suitable mechanism of the insulator- 
conductor transition has been formulated. These results were presented 
during the meeting of the American Physical Society, Dallas, Texas, 

3-12 March. In addition, we are in the process of preparing a manuscript 
to be sent for publication within a few weeks. The findings are summarized 
below. 


The oxidizing dopant removes electrons fi?om a (d^ + dff )-like orbital 

of the insulator, generating an unstable partially oxidized structure. This 
structure is stabilized by a process which results in a decrease of the ^ 
inter-stack distance, accompanied by the energy level reordering of the dz 
ard (d^ + df ) orbitals and partial occupation of the dz2-like orbital. 

Figure 1 shows a contour map of the electron density of the (d^» + df ) orbital 
in the plane of the 3 Pt atoms and 2 CN groups , of the central Pt(CN) 4 2- unit 
of the cluster. Fig. 2 is a similar contour map, but for the dz 2 -like orbital. 

The HOMO in the trimer with the insulator inter-stack separation is 
essentially a non-bonding orbital that is 99.5% localized in the central 
unit. Since there are no evidences of "surface states" as a result of the 
cluster model, this orbital probably represents the true nature of the 
valence band edge in the extended chain. Partial oxidation or removal of 
nq electrons from this band will dreate a heterogeneous distribution of 
charges along the Pt(CN) 4 units. This decrease in negative charges 

of half of the Pt(CN) 4 units must in turn produce a decrease in the inter- 
stack Coulomb repulsion and consequently a shortening of the inter-stack 
distance. Hence, the antibonding dz 2 -like orbital will be pushed up in 
energy, becoming the new valence band edge. At the same time that the 
energy crossing of these orbitals is taking place, homogenization of the 
charge distribution is occuring. Tf*^ excess of positive charge (in relative 
terms) generated every other unit due the removal of electrons from the 
completely localized (d^p + d^» ) MO's is returned to the delocalized dz 2 - 
like orbital. This condition is met by transfer of nq electrons from the 
dz 2 -like to the (d f + d^ )-like band. The resulting valence band is partially 
filled, the charge of each stack being -2 +q. Two competing factors resulting 
from the partial oxidation of the dz 2 -band are principally responsible for 
maintaining the equilibrium prescribed above. Since the dz 2 valence band 
is antibonding between neighboring Pt atoms, its partial oxidation results 
in a strengthening of the Pt-Pt bond. However, further bond-length 
contraction is prevented by the net positive charge q created around each 
Pt atom. 
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Monthly Technical Progress Report 
for the period 

March 23, 1982 to April 22, 1982 

"QV Antmm Chemical calculations for 
Polymers and Organic Compounds" 


A first draft of a manuscript reporting the phenomenon of transition 
from insulator to conductor in tetracyanoplatinate complexes was 
finalized during this period. This paper deals principally with the 
characteristics of the highest occupied molecular orbitgls (HOMO) because 
of their key role in the partial oxidation process. Concurrently we have 
started examining the charge densities originated from MO's located at 
energies below the HOMO. Of special interest is the local bonding along 
the Pt backbone. The charge densities for the trimer (Pt(CN) 4 ] -i - ~ 

Pt-Pt separation of 3.5A have been obtained for several MO'sv Th 
results for three of these MO's are shown in Figures 1 to 3. In all cases 
the contours are expressed in electrons bohrs ' . Their values are: 

0.00001 (0), 0.00005 (1), 0.00025 (2), 0.00100 (3), 0,00200 (4), and 
0.00500 (5). The orbital shown in Figure 1 (with an eigenvalue of 
-0.686 Rydbergs) is d 2 -like, and has bonding character between Pt atoms. 
Figures 2 and 3 correipond to charge densities of (d £ + d'jjJ-like 
orbitals with eigenvalues -0.541 and -0.563 Rydbergs respectively. These 
orbitals have charges localized on the peripheral Pt(CN> 4 units, and 
consequently do not contribute to bonding between stacks. All the other 
MO's studied to date contribute -very little to bonding between Pt(CN> 4 
units. Analyses of the remaining orbitals are ongoing. 

SCF-X ot-DSW calculations leading to the interpretation of photo- 
absorption spectra of aqueous solution of Pt(CN) 4 [1] have been 
initiated. In these calculations we make use of the transition state 
concept introduced by Slater [2] in which the relaxation effects of the 
orbital energies produced by electron excitation are taken into 
consideration. Our calculations are useful in this study since 
relativistic effects have been taken into account fully within the one- 
electron local exchange framework. Comparisons with previous nonrelati- 
vistic work by Interrante and Messmer (3) are being made and will be 
summarized in the next progress, report . 

Consultation with Dr. M. Golub is continuing on quantum chemical 
computations. The principal approach used by Dr. Golub is the semi- 
empirical MNDO molecular orbital method. In conjunction with Dr. Golub, 
discussion of MNDO results on dichloroethyl ether and sulfide have been 
ongoing. 
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